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Some News about Generalised Parton Distributions 

Markus Diehl 
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I briefly discuss some recent developments (and recall some old news) in the theory and 
phenomenology of generalised parton distributions. 
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1 Parametrising generalised parton distributions 




As experimental data on deeply virtual Comp- 
ton scattering (DVCS) and on exclusive meson 
O ■ production is becoming more and more pre- 

cise [2], the demands on an adequate theory 
description are rising. An important part of 
fH ' any theoretical analysis using the framework 

of generalised parton distributions (GPDs) is 

to devise a paramctrisation of these functions. t 
. This is a complex task since GPDs depend on 

(SJ ; two independent momentum fractions x and £ Figure 1: Basic variables describing a GPD. 

t"***" . and on the invariant momentum transfer t (see 

Fig. [J) and since the dependence on these variables is subject to a number of consistency 
constraints. Understanding this dependence is, however, not only a practical necessity but 
can improve our very understanding of how partons are distributed in the nucleon. 

An important constraint on GPDs is the so-called polynomiality property: the nth Mellin 
moment in x is a polynomial in £ with degree n or n — 1 , depending on the particular GPD in 
question. Being a direct consequence of Lorentz invariance, this has turned out to be essential 
even at the practical level, since it is required for the consistency of dispersion relations for 
the processes where GPDs appear (see below). Several methods that ensure this property 

$H ' in the construction of GPD parametrisations have been used in phenomenology: 

' 



• the double distribution representation of GPDs. The Radyushkin-Musatov ansatz J3J [4] 
has been used in a large number of applications — including the VGG code [5] — but it 
should be emphasized that this is only one particular ansatz for the form of double 
distributions. A Polyakov- Weiss D-term [6] is typically added to this ansatz so as to 
provide the term of order £ n in the nth Mellin moment of the distributions H and E. In 
recent work [7] the question was raised whether such a D-term can indeed be freely chosen 
or whether it is already fixed by the double distribution piece due to further consistency 
requirements. 

Recall also that there are several alternative double distribution representations of GPDs 
which generate the term of order £™ in the Mellin moments by themselves and hence do 
not require a D-term from the start [8] . 



*Talk given at the 16th International Workshop on Deep Inelastic Scattering and Related Subjects, 
University College London, England, 7-11 April 2008 [J 
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• the so-called dual parametrisation [9] starts with a partial-wave decomposition in the t- 
channcl and resums these partial waves into a series of two- variable functions (in contrast 
to the double distribution case, where one deals with a single three- variable function). 
In practical applications, only the first or the first two functions in this series have been 
retained. 

Let me note that, in analogy to the magnetic and electric combinations of the Dirac and 
Pauli form factors, the t-channel partial waves are described by the combinations H + E 
and H+j^E of GPDs [10]. This feature is not incorporated in the parametrisations [9]. 

• an ansatz for the Gegenbaucr, or conformal moments of the GPDs [7]. These moments 
depend on £, t and a moment index j, which is to be considered as a variable in the 
complex plane in order to retrieve the GPDs. This method generalises the representation 
of usual parton densities in terms of their Mellin moments, and it has the advantage that 
evolution and a s corrections can be implemented at NLO and partially at NNLO in a 
technically simple and numerically efficient way. 

The lowest Mellin moments of the quark GPDs H 9 and E q are related to the electromag- 
netic Dirac and Paul form factors, which are well measured experimentally and are typically 
used as constraints in GPD parametrisations. Similarly, the axial form factor gives the 
lowest moment of H u — H d . Significant progress is being made in calculating higher Mellin 
moments of quark GPDs on the lattice [TT] . Qualitative features found in these calculations 
provide guidance for parametrising GPDs already now, and eventually one can hope to use 
lattice results as a quantitative input. 

An indirect constraint on the nuclcon-helicity flip distribution E 9 for gluons is provided 
by the sum rule J dx E 9 + J2 q J dx xE q = at t = and £ = 0. Lattice calculations 
give a very small result for J dx xE q , with large cancellations between u and d quark 
contributions. This is in line with the small value of the lowest moment J dx E q inferred 
from the magnetic moments of proton and neutron. Under the assumption that E 9 does not 
change unusually fast when going to finite t and £, the above sum rule hence implies that 
the distribution E 9 is itself small — in contrast to its counterpart H 9 — unless it has one or 
more nodes in x. Current models p21 [13] do not allow for such nodes (for which there is no 
indication from data) and consequently have E 9 much smaller than E u or E d . Going back 
to the point t = and £ = 0, one should note that if / dx E 9 is small then the total angular 
momentum ^ J dx (H 9 + E 9 ) carried by gluons in the nucleon according to Ji's sum rule is 
rather large, because J dx H 9 = J dxxg(x) is about one half. 

Further nontrivial constraints on GPDs are positivity conditions, which ensure the inter- 
pretation of GPDs as densities or interference terms in impact parameter space. As in the 
case of the usual parton densities, positivity conditions are preserved by leading-order evo- 
lution to higher scales but need not hold at NLO or beyond. It seems reasonable to demand 
some explanation in cases where they break down (e.g. in the region of very small x, where 
NLO corrections are known to become important). The most general form of the positivity 
conditions for GPDs is quite involved [14j . Simplified versions concern the skewness effect 
and give e.g. an upper limit on H q (x,£,t) for x > £ in terms of the geometric mean of 
the quark densities evaluated at (x ± £)/(l ± £) Q~51 [3J. To implement or even check these 
constraints in GPD parametrisations is difficult and most often not done. An exception are 
approaches based on calculating double distributions in spectator models [16] . which allow 
both polynomiality and positivity to be implemented automatically, but so far have barely 
been used in phcnomcnological studies. A different type of simplified positivity condition 
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involves the GPDs at £ = and limits for instance E in terms of H and H [T7] . For simple 
analytic forms of the GPDs these conditions can readily be checked — unfortunately this is 
not always done in practice. 

An approximate power-law dependence at small x, as suggested by simple Regge phe- 
nomenology, works well in phcnomcnological parametrisations of the usual parton densities, 
and it is natural to extend this to a behaviour proportional to x~( a+a for GPDs. It is 
important to note that DGLAP evolution changes the value of the effective shrinkage pa- 
rameter a' with the factorization scale [TB], just as it changes the value of the effective power 
a at t = 0. This change can be rather slow, so that the mixing of gluons and sea quarks 
under evolution does not by itself guarantee that a' is similar for these partons at moderate 
scales [TB] . 

2 From exclusive processes to GPDs 

It is well known that, to leading order in a s , the imaginary part of the amplitude for DVCS 
or light meson production involves only GPDs at the special points x = £ and x = — £. 
While the impact parameter representation is somewhat involved for GPDs at nonzero £ 
[19j , it simply gives the distance between the struck parton and the spectator system in this 
case [20] . 

Whereas the factorization formula for the real part of the DVCS or meson production 
amplitude requires the GPDs in the full x-region, a representation based on dispersion 
relations [21] involves, at leading order, only the GPDs at \x\ = £ plus a subtraction term 
which depends on the coefficients of £ n in J dx x n ~ 1 H(x,t;,t) and can hence be expressed 
through the D-term in the Polyakov- Weiss representation. This implies that at tree-level 
accuracy, exclusive processes are only sensitive to GPDs at \x\ = £ and to one t-dependcnt 
subtraction term for each parton species. Both the evolution of the GPDs and explicit NLO 
corrections to the amplitude involve the region |x| > £, in addition to a more complicated 
subtraction constant jTQlEZ] - To access this information experimentally, one must be sensitive 
to violations of Bjorken scaling and thus needs a sufficient lever arm in Q 2 at given xb- 

The polynomiality of Mellin moments, which ties together the GPDs at \x\ < £ and at 
M > £j is crucial for the consistency of the dispersion relations just mentioned. In [7] an 
explicit construction is given that allows one to reconstruct a GPD in the region \x\ < £ 
from its knowledge at \x\ > £, apart from a possible D-term ambiguity. 

Deeply virtual Compton scattering not only offers a large number of observables that 
can be evaluated in the GPD framework [22] but is also under very good theoretical control, 
with radiative corrections known at NLO and in part at NNLO. In numerical studies these 
corrections were found to be of moderate size |23] , except for important effects in scaling 
violation at small x, which are of the same nature as those in inclusive DIS and stem from 
the behavior of the singlet evolution kernels. Nevertheless, NLO corrections of the order of 
20% are not uncommon in collider or fixed-target kinematics. While a leading-order analysis 
of DVCS data should be adequate to reveal basic features and to provide a starting point, I 
see no good reason why analyses should be limited to leading order (including the neglect of 
evolution) as experimental data is becoming more precise. The NLO scattering amplitudes 
for DVCS are reasonably simple for practical use, and the LO evolution of GPDs has been 
implemented in a fast numerical code |24] . 

Meson production is harder to analyse quantitatively, as both power corrections [13] 
and higher orders in a s [25j [12] are larger than for DVCS according to the estimates in 
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the literature. The often substantial size of NLO corrections should not surprise us, given 
that a s ~ 0.4 for scales relevant to a large part of the experimental data. Cross section 
ratios are often assumed to be less sensitive to corrections. The study [12] found that this is 
sometimes true but not always. The ratio of ui and p production cross sections at moderate 
xb is indeed quite stable w.r.t. NLO corrections. For the cross section ratio of and p 
production this is, however not the case, since the former is dominated by gluon exchange 
and the latter by a mixture of quark and gluon exchange, which are affected in different 
ways by NLO corrections. For the GPD models considered in [12], the transverse target spin 
asymmetries for p, cj, and (f> production are very small at tree level, due to cancellations 
between different contributions and to a small relative phase between the amplitudes that 
have to interfere in order to make the spin asymmetry nonzero. The NLO corrections on 
these small asymmetries were found to be large. A careful case-by-case analysis is hence 
required before one can assert the stability of cross section ratios or asymmetries. 
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